Existence of solutions to degenerate parabolic equations via 
the Monge-Kantorovich theory. 



Martial Agueh 

School of Mathematics 
Georgia Institute of Technology 
Atlanta, GA 30332, USA 
agueh@math.gatech.edu 



Abstract 

We obtain solutions of the nonhnear degenerate parabohc equation 
1^ = div{pVc^ [\/{F'{p) + V)]} 

as a steepest descent of an energy with respect to a convex cost functional. The 
method used here is variational. It requires less uniform convexity assumption than 
that imposed by Alt and Luckhaus in their pioneering work [3]. In fact, their 
assumption may fail in our equation. This class of problems includes the Fokker- 
Planck equation, the porous-medium equation, the fast diffusion equation, and the 
parabolic p-Laplacian equation. 
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1 Introduction 



We consider a class of parabolic evolution equations, so-called doubly degenerate parabolic 
equations. These equations arise in many applications in physics and biology [TU] . 
j20j . [22j . They are used to model a variety of physical problems: the evolution of a 
fluid in a certain domain: porous-medium equation ^Sj, Fokker-Planck equation |12j . 
etc. In this work, we focus on these parabolic equations of the form 

' ^ = div (a(6H,Vn)) on {O,oo)xn 

< u{t = 0) = uq on n (1) 

^ a{b{u),Vu) ■ 1/ = on (0,00) x 9 J7, 

where 

a{b{u),Vu) := f{b{u))Vc*[V{u + V)], 
and c* denotes the Legendre transform of a function c : R'^ — > [0,oo), that is, 

c*{z) = sup{(x,z) — c(x)}, 

for z € R'^. Here, is a bounded domain of R"^, v is the outward unit normal to 
Sfi, b : R ^ R is a monotone nondecreasing function, V : Cl ^ R is a potential, 
c : R'^ [0, 00) is a convex function, / is a nonnegative real-valued function, and 
uq : ^} ^ Ris a measurable function. The unknown is u : [0, 00) xQ ^ R, u = u{t, x). 

In a previous work, Alt and Luckhaus [Hj proved existence of weak solutions to (jT]). 
when V = 0, under the following ellipticity condition on a{t,z) := f{t)'Vc*{z), 

{a{t,zi) - a(t,Z2),zi - Z2) > X\zi- Z2 (2) 

for some A > and p> 1, and for all zi, Z2 € R'^. This amounts to imposing that / is 
bounded below, and the cost function c satisfies the ellipticity condition 

(Vc*(zi) - Vc*(z2), zi - Z2) > A I zi - Z2 (3) 

Note that when c{z) = or equivalently c*{z) = where p > 1 is the conjugate of 
g>l, 1 + 1 = 1, condition © reads as 

(I Zi r^Zi - I Z2 r^Z2, Zi - Z2) > A| Zi - Z2 (4) 

which holds only if p > 2. In fact, when 1 < p < 2, the reverse inequality in (jlj) holds 
(see jS], pp. 13). In |3j, the authors approximated by a time discretization, and 
they used a Galerkin type argument to solve the resulting elliptic problems. In the same 
paper, they proved uniqueness of solutions to (0) when V = 0, assuming that ^ holds, 
and the distributional derivative -^p- of a solution -u of is an integrable function. 
The last condition was removed by Otto in ^^1, using the technique of "doubling of 
variables", introduced by Kruzkov [IB], which consists of doubling the time variable 
of two solutions of (0), and treating each solution as a constant with respect to the 
differential equation satisfied by the other solution. 
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In this work, we ehminate assumption (jSJ, and we impose instead, the following 
growth condition on the function c: 

P\z\'i <c{z) <a {\z\'i + 1), (5) 

for z E M'^ and for some a, /3 > and q > 1. Notice that © is much weaker than the 
ellipticity condition ^ imposed by Alt and Luckhaus in 3J. Typical examples are the 
functions c{z) = or equivalently c*{z) = with | + | = 1 and 1 < p < 2. Such 
functions satisfy ((Sj but not © or @ when 1 < p < 2. As mentioned before, they 
actually satisfy the reverse inequality in © or Q when 1 < p < 2. 

We interpret as a dissipative system, and then, we introduce the internal energy 
density function F : [0, oo) R, satisfying F' = b^^. Setting p := b{u), pQ := &(no), 
and /(x) = max(x,0), we rewrite (0) as 

f If + div(p[/p) = on {O,oo)xn 

p{t = 0)=po on n (6) 

pUp-v = Q on (0, cxd)x(90. 

Here, 

Up:=-Vc*[V{F'{p) + V)] 

denotes the vector field describing the average velocity of a fluid evolving with the 
continuity equation in ©, /Oo : — > [0)Co) is the initial mass density of the fluid, and 
the unknown p : [0, oo) x ^ [0,cxd), p = p{t,x), is the mass density of the fluid at 
time t and position x of 0,. The free energy associated with the fluid at time t £ [0, oo), 
is the sum of its internal energy and its potential energy, 



E{p{t)):= / [F{p{t,x))+p{t,x)V{x)]dx. 
Jn 

Problem © includes the 

• Linear Fokker-Planck equation: 

^ = Ap + divipVV) 

I |2 

(c(z) = and F{x) = xlnx) 

• Porous-medium and Fast diffusion equations: 

(V = 0, c(z) = and F{x) = ^ with 1 / m > 1 - i). 

• Generalized heat equation: 

^ = div (\ Vp^ I^^^Vp^ 
(y = 0, c{z) := ^ with | + |, and F{x) = ^x\nx with p > 1). 
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Parabolic p-Laplacian equation: 

(V = 0, c(z) = ^ with I + 1 and F(x) = with m := f^r and p > fi^). 

Doubly degenerate diffusion equation (see [15!): 



^ = div (I vp" r- Vp'^) (7) 

-, and F{x) = with m := n + Sff 



= 0, c(z) := ^ with i + i, and = with m := n + Sff and 

p-1 T 1^ ^ d{p-l) >■ 

The above restrictions on m, n and p are made so that F satisfies the assumptions (HF^ 
and (HFUl) below. 

We are interested in the following questions: under what conditions does © have 
solutions? Is the solution unique? What are the most relevant conditions on c, F and 
y, which ensure that solutions converge asymptotically to an equilibrium? 

In this work, we answered the first and the second questions. We proved existence 
and uniqueness of weak solutions to ©, when the initial mass density pQ is bounded 
below and above, that is, /jo + ^ G L°°(0) (see Theorems 13.111 and 13 . 1 2|) . This restriction 
was made to simplify the proofs, and not to bury fundamental facts into technical 
computations. We include in R,emark l3.131 a method which may be used to extend our 
existence result to the cases where ^ fails to be bounded, and where po belong to a 
wider class of probability densities po £ L^iS^)-, P ^ Q- In a coming paper, we establish 
large time asymptotic results for solutions of ©. 

Our approach in studying existence of solutions to © was inspired by the works 
of Jordan-Kinderlehrer-Otto ^2 and Otto [ISj- In P2]; the authors observed that the 
Fokker-Planck equation can be interpreted as the gradient flow of the entropy functional 



H{p) := [ (plnp + pV) dx, 



with respect to the Wasserstein metric d2- Recall that d2 is a metric on the set of 
probability measures on 12°', with finite second moments, defined by 



inf 



I 1 2 

\x - 



d-f{x,y) : 7 € r(po,Pi)} 



1/2 



where r(pO)Pi) denotes the set of probability measures on R"^ x R'^, having po and pi 
as their marginals (see the definition below). This idea was generalized by Otto in jl5j 
for doubly degenerate diffusion equations of the form 0. 

Let us outline here the proof of our existence theorem to ©• For the sake of 
illustration, we assume that V = 0. The proof consists of four main steps. 
Step 1. We interpret © as a "steepest descent" of the internal energy functional 



Vam3p^Ei{p) := / F(p(x)) dx 

Jn 
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against the Monge-Kantorovich work W^, where /i > is a time-step size, and 'Pa(^) 
denotes the set of probabihty density functions p : 17 — > [0, oo). In other words, given a 
mass density p'^_i of the fluid at time t^-i = {k — l)h, we define the mass density at 
time tk = kh, to be the unique minimizer of the variational problem 

(Pt) : inf {hW^ (pti, s) + E,ip) } (8) 

(see Proposition So, at each time t, the system tends to decrease its internal energy 
Ei(p), while trying to minimize the work to move from state p{t) to state p{t + h). 
Step 2. We write the Euler-Lagrange equation of (-P^), and then, deduce that 



h 



div{p^Vc^[v(F'(p^))]}+Afc(M, (9) 



weakly, for k £ N (Proposition E^J, where Ak{h) tends to 0, as h goes to 0. © shows 

clearly why (jH)) is a discretization of ©• 

Step 3. We define the approximate solution to ®, as 

p^(t,x) = p'^ix) ii te{{k-l)h,kh], ke N 
/(0,x) = po{x), 
and we deduce from @ that, p^ satisfies 

' ^ = div|p'^Vc*[V(F'(p'^))] | + A(/i) on (O,oo)xf] 

(10) 

/(t = 0) = Po on n 

in a weak sense fProDosition l2.9|) . where A{h) is shown to be (/i'^^''^) , e{q) := min(l, g — 
1) (Proposition 

Step 4. We let /i go to in (jTU]). and show that {p^)h converges to a function p, 
which solves (jnj in a weak sense. Here, two convergence results are established: the 
weak convergence of {p^)h to p in ((0, T) x il) , < T < cxd, for a subsequence, 



which proves that y~§rj converges weakly to in ^C^{M x 1?*^)] , and the weak 
convergence of the nonlinear term (divjp'^Vc* [V {F' (p^))]}) to div{pVc* [V {F'{p))]} 
in [C^{R X M'^)]' , for a subsequence. 

The first convergence follows from the second bound in of Proposition 12.31 a con- 
sequence of the maximum principle stated in Proposition 12.21 fsee Lemma l3.3() : starting 
with a probability density function pq which is bounded above, that is, po < N a.e., the 
probability density function p^ - solution of (-P^) - is bounded above, as well, that is, 
p^ < N a.e., for k G IV. As a consequence, {p^)h is bounded in ((0,oo) x 17), and 
then, converges to some p in ((0, T) x oo) , < T < oo, for a subsequence. 
The second convergence is one of the most difficult tasks in the proof of the existence 
theorem. Its proof requires elaborated intermediate results. Here, we see some techni- 
cal differences with the works in 3.^ and |15j . Indeed, due to the weaker condition © 
imposed on c here, compared to the stronger ellipticity condition @ or @ in jSI, the 
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method used in j2| and do not yield, here, strong convergence of the nonhnear term, 
as in IS] and ^H]- So, to prove the - weak - convergence of the nonhnear term, here, we 
proceed as foUows: 

(i) . First, we improve the previous convergence, by showing that, in fact, {p^)h con- 

verges strongly to p, for a subsequence, in ((0,T) x $7) (Proposition I.S.7|) . 

(ii) . Then, we deduce that (div{/Vc* [V converges weakly to 

div{pVc* [V{F'{p))]} in [C^{R x R'^)]' , for a subsequence fTheorem mill) . 

To prove (i), one needs to have a good control on the spatial derivative of p^, for example, 

to show that | V | is bounded in L^* ((0,r) x Q) , < T < oo. The main 

ingredient used to establish this result is the following Monge-Kantorovich type energy 
inequality: 

E,{p~o) - Ei{p\) > / (V , ~S{y) - y)p~i{y) dy, (11) 



for pq, jji £ Va{^)- Here, S denotes the c-optimal map that pushes pi forward to po (see 
the definition in Proposition II. 1|) . A more general statement of the energy inequality is 
given in Theorem 12.81 (fTTj) can be seen as a consequence of the displacement convexity 
of the internal energy functional 'Pa(^) 3 P^ Ei{p)^ that is, the convexity of 



[0,1] Bt^E,{pi_t) 



where. 



~pi-t ■■= (1 -t)id +tS 



# 



Pi 



(12) 



is the shortest path joining pi and po in Va{^)- When c{z) = in which case S is 
the gradient of a convex function, the interpolation in (|12|) was introduced by McCann 
in [T?| . 

Indeed, setting po := p^_i and pi := p^ in (fTT|). and using the Euler-Lagrange 
equation of (-P^ ), that is. 



h 



^[F'ipt: 



(13) 



where 5^ is the c (^) -optimal map that pushes p^ forward to p^_i, we obtain that 

h I {v{f\pI)),Vc*\v{f\pI))])pI<E,{pI_,)-E,{pI). (14) 



We integrate (|14|) over t G [0, T], and we use Jensen's inequality, to deduce that 



Jn 



(V(F'(/)),Vc'^ V(F'(/)) )p^<EM-\^\F 



(15) 



We use condition © - precisely, c{z) < a {\ z + 1) -, combined with (fT5|) and the fact 
that (p'^)/i is bounded in L°° ((0, oo) x $7), to conclude that 



T 



Jn 



VlF'ip 



< est (see Lemma ESI)- 
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Then, we use that is bounded in L°° ((0,cx3) x Q,) (see (|26|) ) - a consequence of 

the minimum principle of Proposition 12.21 ~ to deduce that |v (-F'(/o'^)) | is bounded 
in L-?* ((0,r) X f)) , < T < oo. This yields (i). 

To prove (ii), we first use (|T3|) and condition © - precisely c{z) > P\z\'^ - to have 
that |Vc* [V(F'(/))] } is bounded in L'^ (0 x (0,oo)) fLemma 13. 8() . from which we 

deduce that |Vc* [V | converges weakly to some a in L"^ (il x (0, T)), for 

a subsequence, and for all < T < oo. Next, we use (i) and the boundedness of 
|v {F'{p^)) 1^ in L?* {fl X (0,oo)), to obtain that |v |^ converges weakly to 

V {F'(p)) in L"? (il X (0, T)) for a subsequence (Lemma 13.8(1 . In the end, we extend the 
energy inequality (|1H) in time-space (Lemma 13. 9|) . and we combine the new inequality 
with the strong convergence of {p'^)h to p, the weak convergence of |v (^F'{p^)^ | to 

V(-F'(p)), and the weak convergence of {VC^ [V }^ to fj, to establish that 

(divjp'^Vc* [V })^ converges weakly to div{pa) for a subsequence, and that 

div(/9<T) = divipVc'' [V (F'(p))]} (Theorem l3Tnjl . The convexity of c* plays an impor- 
tant role in this proof. 

Notations 

• n c R'^, d>l is open, bounded, convex and smooth, and ^It ■= iO,T) x ft, for 

< T < oo. 

• Bji{x) C R'^ denotes the open ball of radius R, centered at x, Bji{xY '■= R'^ \ 
Bji{x), and p* denotes the conjugate index of p > 0, that is, | + ^ = 1. 

• Va{^) := |p : ^ [0,cxd) measurable, J^p{x)dx = l|, and for < i? < 
oo, Vi^\n) := |p G Va{n) : p<R a.e. }. 

• If If : Q ^ R, then || <yJ ||_L9(n) denotes the L'^-norm of if, and spt (99) denotes the 
support of If, that is, the closure of {x £ Q : ip{x) ^ 0}. 

• If 2; = (xi, • • • , Xd) a nd y = (yi, • • • , ^d) are vectors in R'^, then (x, y) := Ya=i ^iVi^ 
and I X I := ^J {x, x). 

• If A C R^, D > 1, is convex, and G : A ^ 1? is convex, then G* : R^ ^ R 
denotes the Legendre transform of G, that is, 

G'(y) := ^up {(X, ,) - <}(.)}, where G(.) := { °W « >1 

By abuse of notations, we will identify G and G. 

• If j4 C R'^ is Borel, then | A \ denotes the Lebesgue measure of A, and Ia{x) '■= 

f 1 ii X G A denotes the characteristic function of A. 

1 U otherwise, 
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Throughout this manuscript, M and N are positive reals, a.e. refers to the d-dimensional 
Lebesgue measure, and 



Definitions 

Probability measures with marginals. Let /Uq and fii be probability measures on 
M"^. A Borel probability measure 7 on the product space M'^ x R'^ is said to have //q 
and Hi as its marginals, if one of the following equivalent conditions holds: 

(i). for Borel A C R"^, 



(ii). For {ip, il)) G L^^ {R'^) x L^^ (-K*^), where L^. {R'^) denotes the space of /Lij-integrable 
functions on R'^ {i = 1,2), 



We denote by T(ij.q, fii), the set of all probability measures satisfying (i) or (ii). If fiQ 
and /Il are absolutely continuous with respect to Lebesgue, and po, pi denote their re- 
spective density functions, we simply write r{pQ,pi). 

Push-forward mapping. Let pq and pi be probability measures on R'^. A Borel map 
T : R'^ —>■ R'^ is said to push po forward to pi, if 

(i) . pi[A]= po[T~^{A)] for Borel A C R'^, or equivalently 

(ii) . ip{y)dpiiy) = lR^V>iT{x)) dpoix) for p G 1.^,(1?^). 

Whenever (i) or (ii) holds, we write that pi = T^po, and we say that T pushes pq 
forward to pi. 

The next proposition is due to Caffarelli |5, and Gangbo-McCann IQ'. It asserts the 
existence and uniqueness of the minimizer for the Monge-Kantorovich problem. 

Proposition 1.1 (Existence of optimal maps). 

Let c : R'^ — > [0, 00) be strictly convex, and po,pi G Va{^). Then, 

(i). there is a function v : 0, ^ R such that, T : Cl ^ Cl, T := id — (Vc*) o Vn pushes 



po forward to pi, where u{x) = miy^-^ |c(x — y) — v{y) | for x G 
(ii). T is the unique minimizer (a.e. with respect to po) of the Monge problem 




and 




-f[A X R'^] = po[A] and j[R'^ x A] = pi[A]. 
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(in). The joint measure 7 := {id x T)^pQ uniquely solves the Kantorovich problem 
(/C) : inf I / c{x - y) d'y{x,y), £ T{po, pM . 

(iv). T is one-to-one, that is, there exits a map S : 0. ^ pushing pi forward to po, 
such that T {S{y)) = y a.e. with respect to pi, while S (T{x)) = x a.e. with respect 
to pq. 

Moreover, S = id -\- Vc*(— Vu), where v{y) = inf^g^ |c(x — y) — u(x)| for y G $7. 
V is called the c-transform of u, and it is denoted by v := u'^. 

We will refer to T (respectively S) as the c-optimal map that pushes po (respec- 
tively pi) forward to pi (respectively po), and 7 will be called the c-optimal measure in 
r(/3o,Pi)- 

Wasserstein metric. Let c : -R*^ — > [0,oo) be strictly convex, h > 0, and po, pi G 
Va{^). We define 

W^{po,pi) := inf {y^^^^^c(^) d7(x,y): jeT{po,pi)}. 

If c{z) = we denote Wj} by . When c{z) = ^ and /i = 1, da := ^Jw^ is called 
the Wasserstein metric. 

We deduce from Proposition II. II that, there exist a unique probability measure 7 G 
r(po5 Pi), and a unique mapping T that pushes po forward to pi, whose inverse S pushes 
pi forward to po, such that 

W,'{po,Pi)= I d7(x,y) = / Po(x)dx 



c 



Assumptions 

(HCl) : c : J?"^ ^ [0, 00) is such that = c(0) < c{z) for z ^ 0. 
(HC2) : lim|^.|^oc = 00, i.e. c is coercive. 

(HC3) : Plzl" < c{z) < a (I z |9 + 1) for z G R'^, where a, /3 > and g > 1. 

(HFl) : Either lima;^+oo — ^ = +00, (i.e. F : [0, 00) M has a super-linear growth 
at +00), or lim^^_|_oo = and F'{x) < 0, Vx > 0. 
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(HF2) : (0, oo) 3 x >—>■ x'^F{x is convex. 

We impose assumption (HF^ to ensure that the Legendre transform F* of F is finite on 
(-00,0). In fact, (HFlU), the strict convexity of F, and F(0) = 0, imply that F*{x) G R 
for X < 0. Then, F* {F'{a)) is finite for a > 0. 

The following cost and energy density functions satisfy the above assumptions: 

• c{z) = "^^^i Ai \ z\'^^, where n G iV, qi > 1, and Ai > 0, (take q = max|j=i ... „}.(gj) - 
qio for some G {1 • • • , n}, l3 = Ai^, and a = Y^'^^i Ai ). 

• F{x) = xlnx, F{x) = m>lorl-i<m<l, and F{x) = ^.7=1 ^ii^) ^ 
where n £ N, A^ > 0, and the Fi are like the previous F. 

2 Calculus of Variations on Va{^) 

We discretize (jSJ, and prove in section 2.1 that 

(P): inf {l{p):=hW^{po,p) + E,{p): p G Vain)] (16) 

admits a unique minimizer pi. The reason why we minimize such a functional I{p) will 
be clear in section 2.2, where we find the Euler-Lagrange equation of (P). In fact, we 
shall see that the Euler-Lagrange equation is nothing but the discretization of In 
section 2.3, we show that 

EM - EM > ^^^l^^, (17) 

where pi-t (|12jl denotes the probability density obtained by interpolating pQ and pi 
along the "geodesic" joining them in Vai^,). We refer to ()17() as the (internal) energy 
inequality. We shall see later on, that (\17\i is an essential ingredient in the proof of 
the convergence of the approximate sequence {p^)h (see the definition in section 2.4) to 
solutions of ©. 

2.1 Existence of solutions to a minimization problem (P) 

Throughout this section, h > 0, and pq G 'Pa(^) is such that po < M a.e. We show that 

(Pr): ini {l{p):=hW^^(po,p)+E,{p): pGP(«)(l^)} (18) 

admits a unique minimizer piR, for R>M (Proposition EU, and that pm G vi^"^\n) 
for R > 2M, that is, < pm < M a.e. (Proposition E^l) . We deduce that (P) ^ has 
a unique minimizer pi, which satisfies < pi < M a.e. (Proposition . 

Proposition 2.1 Let R > M, and assume that F : [0, cxd) R and c : R'^ —>■ [0, 00) 

are strictly convex, and c satisfies (Hd^. Then (Pr) has a unique minimizer piR, which 
satisfies 

\n\F < E,{pir) < Eiipo). (19) 
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Proof : Let linf denote the infimum of I{p) over p € Va {^)- Since pQ € Va {^), 
Ei{po) < oo and c(0) = 0, we have that 1^/ < j^Ei{po). Moreover, because of Jensen's 
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inequaUty and the fact that c > and p G Va{^), we have that linf ^ yjIT 
We deduce that linf is finite. Now, let (p^^^)^ be a minimizing sequence for (Pr). We 
have that is bounded in L°°(Q). As a consequence, {p^"'^)n converges weakly-* 

to a function piR in L°°(J7), and then, weakly in for a subsequence, since O is 

( R) 

bounded. Clearly, piR € Va {^)- Furthermore, because of Proposition 5.3.1 2', we 
have that Va{^) B p ^ I{p) is weakly lower semi-continuous on L^{Q), as the sum of 
weakly lower semi-continuous functions. Therefore, 

I{Pir) < liminf = linf < I{pir), 

n—foo 

which shows that piR is a minimizer of (Pr)- The uniqueness of piR follows from the 
convexity of Va{^) 3 p ^ Wl!;{po, p) and the strict-convexity of Va{^) ^ P ^ P{p) 
(see Proposition 5.3.1 0). 

Next, we observe that I{pir) < I{po), and since Wl}{po, po) = and W^^(po) Pii?) ^ 
(because of (HCQ), we deduce that Ei{piR) < Ei{po). We use Jensen's inequality and 

the fact that piR £ V^\^), to conclude that | | F (^yjfi-^ < Eiipm) D 
Proposition 2.2 (Maximum/minimum principle) 

Let R > 2M , and pQ be such that N < po < M a.e. Assume that F : [0, oo) M and 
c : — > [0, oo) are strictly convex, and c satisfies (H(^^. Then, the minimizer piR for 
(Pr) satisfies N < piR < M a.e. Therefore, piR does not depend on R. 

Proof: The proof we present here is similar to that in ^7], where c{z) = and 
F{x) = xln(x). Since the proof of "/Oi-R ^ ^ a-e-" is analogue to that of "pi/J ^ M 
a.e.", we only prove that piR < M a.e. Suppose by contradiction that E := {y £ 
^ '■ PiRiu) > ^} has a positive Lebesgue measure. The idea is to come up with 
Pi^R ^ 'Pi^\^), such that I{pir) > /(/j^'^). This contradicts the fact that piR is the 
minimizer of / over Va^\o,). 

Let "jR be the c/j-optimal measure in T{pq,pir). We have that 

^RiE'' xE)>0, (20) 

where E'^ := \ E; otherwise 

M\E\< [ piR{y) dy = jRiR^ x E) = ^r{E x E) < jr{E x R'^) 
Je 

= / po{x)dx <M\E\, 
Je 

which yields a contradiction. Consider the measure v := ^rIe^xE defined by 



i{x,y)dv{x,y) = I ^(x, y) d7ij(x, y), 

R'^xR'' Je^xE 

for ^ G Cq{R'^ X R'^), or equivalently 

z.(F)=7fi[Fn(^^xS)], 
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for Borel sets F C M!^ x R^. Denote by vq and ui its marginals, that is, 

/ Wix) +'ilj{y)] diy{x,y) = ip{x)duo{x)+ i}{y) dui{y), 

for if, tp & Cq{K''). Since v « and G r(po,/OiR), we have that uq « pQ{x)dx 
and v\ « piR{y)dy. As a consequence, uq and ui are absolutely continuous with 
respect to Lebesgue. Denote by vq and vi their respective density functions. We have 
that, 

(i) . < Vq < M a.e., and < vi < R a.e., and 

(ii) . Vq = a.e. on E, and t^i = a.e. on E"^. 



For e G (0, 1), we define := piR + t{vo — vi), and the probability measure ■j)^' by 
/ ^{x,y)d'y^^\x,y) := / ^(x, y) d7i?(a;, y) 

JRdxRd JRdxRd 

+ e / [^(x, x) - ^{x, y)] d^R{x, y), 

JE'^xE 

for ^ G Co(l?'' X R^). Because of (i), (ii) and the fact that 2M < R, we have that 
< p^il < R, and 

/ p^^liy) dy = 1 + e ^E' x E) - 7«(S^ x E)] = 1. 
Jn 

Hence, G 'Pi^\o,). Moreover, since jr G r(po,PiR), and u has marginals vq = 
vo{x) dx and vi = vi{y)dy, we have that 7^^^ G r(po,Pi^)- Now, we show that I{Pir} < 
I{pir), for e small enough. Indeed, 



(^) 



7(pg) - /(pir) = h te(po, pS) - W^{p^, p^r)] + / [F(pg) - F{p^r) 
Because 7^^-' G r(po, p^^) and c(0) = 0, we have that 



(21) 



wKpo,P'fR)-W^{po,piR) < [ 

JR 



R'^xR''^ 



—f^ ) dj)i'{x,y) 



-L 



Rd-xRd 



el c 

'e'^xE 



x-y 
h 

x-y 
h 



dlRix,y) 
d'yR{x,y). 



On the other hand, according to (i) and (ii), we have, for e small enough, that 



PiR = PiR -evi> M -evi>0 on E, 



and 



PiR = PiR + €Vo >€vo>0 on ET][vo>0]. 



(22) 

(23) 
(24) 
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We combine (i), (ii), (|^ . (|^ . and the fact that F G C^((0, oo)) is convex, and 
u = 'JrIe'^xE has marginals uq = vo{x) dx and ui = vi{y) dy, to obtain that 



n 



< e 



< e 



[ F {piR + evo)- F{pir) ]+ [F {piR -evi)- F{pir) ] 

Je 

F'{piR + €Vq)vq- / F' {piR - evi)vi 
Je 



£;m[tio>o] 



/ F'{M + evo)vo - / F\M - evi)vi 
Je<^ Je 

F'{M + evo{x)) - F'{M - evi{y))) djR{x,y) 



E'^xE 



And since F £ C'^ ((0, oo)), (i) and the above estimate give that 



/ \Fip[%-Fip^R)] =0(6 
Jn L J 



(25) 



Combining (|21l) . (|22j) and H25j) . we conclude that, for e small enough, 



E<^xE 



x-y 
h 



d'yRix^y) < 0, 



where, the last inequality holds because of (HCQ and ()20|1 . 



□ 



Proposition 2.3 (Existence and uniqueness for (P)) 

Assume that N < pQ < M a.e., and F : [0,oo) — > R and c : [0, cxd) are strictly 

convex, and c satisfies (H(^^. Then, pi := piR (defined in Proposition \2.2]) is the unique 
minimizer for (P) ilb}) . Therefore, 



and 



n\F 



N < pi<M a.e., 
1 



n\ 



<E,{pi)<E,{po). 



(26) 



(27) 



Proof : Let p G Va{^)-, and i^P^^^) ^y2M ^ sequence in Va {^) converging to p, such 
that, 



(28) 



as in Proposition 1.4.1 [2]. Since pi is the minimizer for {Pr) f Proposition 12 . 2|) . we have, 
using ((2H1), that 



hW^{p^,pi)+ [ Fipi)<hW^{po,p^''^)+ [ F{p). 
Jn Jn 



(29) 



And, since (p^^^) ^ converges to p in L^{^}), Proposition 5.3.2 [2] gives that 



]imWNpo,p^''A=W^{po,p). 

R]oo V / 



(30) 



13 



We let go to oo in (|29jl . and we use (|3Ujl . to conclude that pi is a minimizer for (P). 
The uniqueness of the minimizer follows from the strict convexity of Va{^) B p ^ I{p) 
as in Proposition 12.11 and (|26|) and (|27jl are direct consequences of Proposition 12.21 and 

(m □ 

In the remaining of this section, we state two propositions needed to establish the con- 
vergence of the approximate solution p^ ()65() of ©, as /i goes to 0. The first proposition 
stated below, shows that the interpolant densities pi-t, t E [0, 1] (fT2|) between two 
probability densities po and pi, which are bounded above, are also bounded above. 

Proposition 2.4 Let po, pi £ Vai^) be such that po, pi < M a.e., and assume that 
c : [0, oo) is strictly convex, of class C^, and satisfies c(0) = and (H(J^. Denote 

by S the c-optimal map that pushes pi forward to pQ, and define the interpolant map 

St := (1 -t) id+tS, 

/or t e [0, 1]. Then, for i € Cc{R^), ^ > 0, 



i {St{y)) pi{y)dy < M / ^x) dx. (31) 
n JRd 

Proof: The proof will be done in two steps. In step 1, we prove (|31|1 for sufficiently 
regular cost functions. Here, we use the fact that V5 is diagonalizable with positive 
eigenvalues when c, c* G C^{M'^), and A (det)^^'^ is concave on the set of d x d 
diagonalizable matrices with positive eigenvalues. In step 2, we approximate a general 
cost function c by regular cost functions c^, and we obtain ()31() in the limit as k goes to 
oo. 

Step 1. c is strictly convex, and c, c* G C^{R'^). 

Proposition 14.11 gives that pi-t '■= {St)#pi is absolutely continuous with respect to 
Lebesgue, for t G [0,1]. Let pi-t denote the density function of pi^t- Then, H31() reads 
as 



i{x) pi^t{x)dx < M I ^{x)dx. 
n Jr<i 

Thus, it suffices to show that pi^t ^ • 

Recall that, because of Proposition 14. H there exists a set K d oi full measure for 
A*i := Pi{v) dy, such that St is injective on and iov y £ K and t G [0, 1], VS{y) is 
diagonalizable with positive eigenvalues, and 

/ Pi{y) = Pi-t (Stiy)) det [VStiy)], (32) 

where VSt{y) = (1 — t)id + tVS'(y). Since po,pi < M a.e., and S^pi = po, we can 
choose K, such that pi{y), po {S{y)) < M foi y G K. We use t = 1 in (|32p. and the fact 
that Po (S(y)) < M, to deduce that 

det[V5(y)]>^. (33) 

Because A i— > (detA)^^^ is concave on the set of d x d diagonalizable matrices with 
positives eigenvalues, we have that 

[ det VSt{y)f' > (1 - t) + t ( det [VS{y)] f'' . (34) 
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We use (|33j) . (j34p . and the fact that pi{y) < M, to obtain that 

det[VSt{y)]>^. (35) 

We combining and we use the fact that St is injective on K, to deduce 

that pi-t < M on St{K). But, since pi{K'^) = 0, and /ii_t = {Stj^pi, we have that 
/ii_t [(5t(K))'^] = 0, and then pi_t = on [St{K)f. We conclude that < M. 
Step 2. c satisfies the assumptions of the proposition. 
Let {ck)k be a sequence of strictly convex cost functions satisfying 

Ck, 4 G CHR''), 

Ck ^ c locally in C^{R'^), as A; ^ cx), (36) 
= Cfc(O) < Ck{z) for z^O 

(see Proposition 1.3.1). Denote by Sk the Cfe-optimal map that pushes pi forward to 
po, and set 

:= (l-t)id + t5,., 

for t E [0,1]. Lemma 2.2.2 [2] gives that ('S'l*^^^ converges to St a.e. on [pi ^ 0] for a 
subsequence, and because of step 1, we have that 

/ ? (S^\y)) Piiv) dy<M [ ^(x) dx. (37) 

We let A; go to oo in (|H7|). and we use that < ^ G Cc{R'^), and Fatou's lemma, to 
conclude (|3T|) □ 

Next, we state a proposition needed in the next section, to prove the strong con- 
vergence of the approximate solutions {p^)h of © in ((0, T) x fi), for < T < oo. 

Proposition 2.5 Lei f : R ^ R and g : R ^ R be strictly convex, of class C^{R), 
such that lim^^oo = oo. Given M, S > 0, and p, q > 1, define 
Am,s ■■= |(^ti,M2) e i^(^2)^ : \\uj \\Li{n) < M, \\g'{uj) \\w^'P{n) < M, and 
In [/'M - nm)][u2 -ui]<5, (j = 1,2)}, 

and set 

Am(5) := sup II -"2 - lii 

Then 

limAM('5) = 0. 

Proof: Suppose by contradiction that there exist k > and {uj)sio-, (j = 1,2), such 
that {uf,U2) G Am,5, and 

II 4 - ^1 llLi(f^) > ('^^) 

By the Sobolev embedding theorem, (^g'{uj)^ converges strongly in LP(J7), and then, 

a.e., for a (non-relabeled) subsequence. Since g G C^{R) is strictly convex and has a 
super-linear growth at oo, we have that {g')~^ is continuous. We deduce that 
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(i). {u^i)s converges to some function uj a.e., for j = 1,2. 



We use (i), || Uj ||l<j(q) < M, and the fact that g > 1 to conclude that {u^j)s converges 



strongly to uj in L^{i}). And since || — M2 llLi(f2) > obtain that 



> K. 



(39) 



Now, we use (i), the convexity of /, and the fact that [ f'{u2) — f'{u\) ] [ "^2 ~ ^1 ] — 
to have that 



< / [nu2) - /'(ni)] [U2 - ui] < hininf / f'{4) - f{u{] 

This implies that 

[/' (^^2(2;)) — /' (ni(x))] [u2{x) — ^1(2;)] = for a.e. x G fi. 



< 0. 



(40) 



Since / G C^{R) is strictly convex, we have that /' is one-to-one, and then, ()4U() implies 
that ui{x) = U2{x) for a.e. x G f]. This yields a contradiction to H39() □ 

2.2 Properties of the minimizer for (P) 

We establish the Euler-Lagrange equation for (P) H16|) . and we derive some properties of 
the minimizer for this problem. The next proposition is the first step towards showing 
that (P) is a discretization of or in other words, © is the steepest descent of the 
internal energy functional Pj, with respect to the Monge-Kantorovich work W^. 

Proposition 2.6 Let po G Pa(^) be such that N < pQ < M a.e. Assume that F : 
[0,00) R is strictly convex, and satisfies F ^ C"^ ((0, 00)), and c : M'^ [0,oo) is 
strictly convex, of class , and satisfies (H(]^ - (H(J^. If pi denotes the minimizer 
for (P), then the followings hold: 



(Vc 



nxn 



x-y 
h 



i'iy)) djix,y) + / P{pi{y))div'ip{y) dy = 0, 



(41) 



for ip G C^(0,1?'^); here P{x) := Pf(x) := xP'(x) - P(x) for x G (0,oo), and 7 is the 
Ch-optimal measure in T{po,pi). Moreover, 

(i). P(pi) G H^i'°°(17). 



(ii). If S is the c^-optimal map that pushes pi forward to po, then 

Siy) - y 



h 



Vc^ [V(P'(/,i(y)))] , 



(42) 



for a.e. y €z Q, and 

pi{y)- pojy) 

h 



ipiy)dy + pi{y){Vc^[V{F'{pi{y)))],V^{y))dy 



< 



2h 



sup 



D'^ip{x) 



x-y|^d7(x,y), (43) 



nxn 



foripe c^in). 



16 



Proof: Since c G C^{R'^) is strictly convex and satisfies (HQ2I), we have that c* G 
C^{R'^), and (Vc)~^ = Vc*. Fohowing we consider the flow map {cl)e)eeR in 

C°°{n,n), defined by 

(po = id. 



(44) 



where ip E Cc°°(0, IZ"'). We have that det (V0e) / 0, and 



d{det V(/>e 



e=0 



div0. 



(45) 



We define on Q, the probability measure /x^ := {4>e)#Pi- Since c/i^ is a C^-diffeomorphism, 
then is absolutely continuous with respect to Lebesgue. Let pe denote its density 
function. Clearly, G Vai^)-, and 



(pe ° 0e) det (V(/)e) = pi a.e. 
Next, we define on x Q, the probability measure 7^ := {id x (/)e)#7, i.e. 



(46) 



i{x,y) d7,(j;,y) 



e(x,0,(y)) d7(x,2/), VeGC(OxO). 



We have that 7^ G r(pO)/Oe), and then, the mean- value theorem gives that 

W^{pO,Pe)-W^{pO,Pl) 



< -[ch{x - My)) - Ch{x - y)] d'y{x,y) 



X - y + 9{y - My))] 



-(y) ) d7(x,y), 



where 9 G [0, 1]. Because of H44() . we have that | | < || V' for e > 0. Then, we 

use that c G C^{R'^), the Lebesgue dominated convergence theorem, and (|lHl . to obtain 
that 

W^i'(PO,P.)-^e^(PO,Pl) , 

iim sup < - 



On the other hand, because of (|l6|l . we have that 



Fip,{x))dx= / F(/5,o0,(y))detV./>e(?/)dy 



(Vc,,(x-y),'0(y))d7(2;,y). (47) 

pi{y) 



^' detV0,(y) 

And since, F £ ((0, 00)), we deduce by the mean- value theorem that 
Fip,{x)) - F{piix)) 



detV(/),(y) dy. 



■ dx 



1 



Pi 



det V</)f 



+ 



Pi 



F(pi) detVcpe + F{pi){detV(l)e - 1) 



det V(/>, - 1 



det V(/)e - 1 



Pi Pi 



(48) 
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where 9 G [0, 1]. We combine (011), (H^J and (gSl) to have that 



hm 

We use (jl7|) and (gHl), to conclude that 



■dy 



P{pi{y))divip{y)dy. 



(49) 



(Vc/,(x-y),V(y))d7(x,y) + i / P(pi(y)) div dy < 0. (50) 

Since Vc/j(z) = Vc(|), and tp is arbitrarily chosen in C^{i},R'^), (jSUl) implies (HTJ. 
(i). By N < pi < M a.e., and since F G ((0, oo)), we have that P{pi) G L°°(l^). 
Now, let (/9 E C;?^(f^), and for an arbitrary i e N, define -0 = (?/'j)j=i,...,d G C^{n,R'^) 
by V'j <5ij V) where (^jj denotes the Kronecker symbol. Because of (j41|) . we have that 



P(pi(y))|^(y) 



dc ( X — y 



< sup 



QxQ dzi \ h 
dc ( X — y 



dzi \ h 



< M W'^Wh^in) sup 



(^(y)d7(x,y) 

95(y) |pi(y)d2/ 
9c ( X — y^ 



dzj V 



And since c € C^{R'^), we deduce (i). 

(ii). Because P{pi) G VF"'^'°°(r2), we can integrate by parts in (|1T|). We use that 7 G 
r(/9o>Pi) and S^pi = po> to obtain that 



{V[P{pi{y))],^{y))dy 
Pi(y)(V[F'(pi(y))] ,V'(y))dy, 



for i) E C7c°°(J^,l?^). And since ■0 is arbitrarily chosen, we deduce that 
'^4^^P^) Pi(y)=V[F'(pi(y))] pi(y). 



(51) 



for a.e. y G 0. We combine (|^. and the fact that (Vc)^^ = Vc* and pi 7^ a.e., to 
conclude (|i2|l . 

Next, we consider (/? G C^(il), we take the scalar product of both sides of with 
pi{y)V(p{y), and we use that 7 = (id x S)^ pi, to obtain that 



nxn 



{y - x,\7(p{y)) dj{x,y) 



( Vc'^ [ V {F' (pi(y)) + Viy)) ] , V <^(y) ) pi (y) dy. 

(52) 

Now, we express ^ Inxniy ~ ^' '^viv) ) djix, y) in terms of /^^ pi iv) ^Poiy) ^p(^y'j ^y. Since 
7 G r(pO)Pi)) we have that 

pi{y) - Po{y) 



h 



nxn 



iviy) - d7(x,y). 
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Combining the above equality with the first order Taylor expansion of 99 around y, we 
obtain that 

t/ {y - x,Vf{y))d-f{x,y) - ]- [ ipi{y) - Poiv)) fiv) 
Jnxn " Jn 

< 777- sup I L>V(a;) I / \x - yf d'y{x,y). (53) 
'^'T' xen Jnxn 

We substitute ^ into ^ to conclude (031) □ 
2.3 Energy inequality 

We establish an inequality relating the internal energy Ei[pQ) and Ei[pi) of two proba- 
bility density functions pQ and pi. This inequality will be called energy inequality and 
will be used later on, to improve compactness properties of the approximate sequence 
p^ (see the definition in section 2.4), to solutions of ©• First, we prove this inequality 
for smooth cost functions c and c*, whose Legendre transform c* are C^. Instead of 
using the density function F, we consider a more general function G, which satisfies 
some assumptions to be specified later on. The (internal) energy inequality reads as 

/ G(po(y))dy- / G{pi{y))dy>- [ Pcipiiy)) div (Siy) - y) dy, (54) 
Jn Jn Jn 

where S is the c-optimal map that pushes pi forward to po, and Pcix) := xG'{x) — G{x). 
For smooth cost functions c, this inequality is simply a consequence of the displace- 
ment convexity of Vai^) ^ p ^ J^G{p{x)) dx, that is, the convexity of [0,1] 9 t 1— > 
G{pi-t{x)) dx, where is the probability density obtained by interpolating po and 
pi along the "geodesic" joining them in Vai^) (see Proposition 14. 1|) . To prove (|54|). 
we rather follow a more direct procedure, using the following result of Cordero and 
Otto (Proposition 14. Ij) : if po, pi G Va{^), c, c* G C^{R'^), and S is the c-optimal map 
that pushes pi forward to po, then VS{y) is diagonahzable with positive eigenvalues for 
Pi '■= Pi{y) dy - a.e. y £ Q. More! over, the pointwise Jacobian detVS, satisRes 

0^pi{y) = detVS{y)po{S{y)), (55) 

for pi - a.e. y £ Q. 

Proposition 2.7 (Energy inequality for regular cost functions) 

Let Po, pi £ Va{^) be density functions of two Borel probability measures po and pi 
on M^, respectively. Let c : R'^ — > [0, cxd) be strictly convex, such that c,c*£ C^{R'^). 
Let G : [0, 00) R be differentiable on (0, 00), such that G(0) = 0, and (0,cxd) 3 x^ 
x'^G{x~'^) be convex and nonincreasing. Then, the internal energy inequality \54\ ) holds. 
Ln addition, if Pg{pi) £ W^'°°{^) and pi > a.e., then 

[ Gipo{y))dy- [ G{pi{y))dy> [ {V[G' {pi{y))\, S{y) - y) pi{y) dy. (56) 
Jn Jn Jn 

Proof : Set 

A{x) := x'^G{x-'^), X £ (0,cx)). 
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We observe that 

A'{x) = -dx^^-^Pcix-^). (57) 
Since A is nonincreasing, we have that Pg > 0, and then 
(i). (0, cxd) B X ^ ^^j^ is nondecr easing. 

Proposition 14.11 gives that VS'(y) is diagonahzable with positive eigenvalues, and that 
holds for /ii - a.e. y G $7. So, po{S{y)) / for pi - a.e. y G $7. We use that 
G(0) = 0, = po, and (|55|) . to deduce that 



G{po{x)) dx 



G{po{x)) 
po^o] Po{x) 

pi{y) 



poix) dx 



G{po{S{y))) 



G 



detVS{y) 



n PoiSiv)) 
detVS{y) dy. 



Piiv) dy 



Comparing the geometric mean (det \/S{y))^^'^ to the arithmetic mean 



trygfa) 



(58) 



we have 



that 



pi{y) 



detV S{y) 



> pi{y) 



d 



trVS(y) 



Then, we deduce from (i) and the above inequality, that 



A 



pi(y)Vrf 



(59) 



where. 



A 



trVSjy) 
d 



Now, we use ()57() and the convexity of A, to obtain that 



pi{y)A 



A 



pi(y)Vrf 



> pi{y) 



pi{y) 



A 



pi(y)Vrf 

Gi {pi{y)) 



+ A! 
d{K-l 



A- 1 

PiiyY''') \pM^^ 
PciPiiy)) 



pi{y) pi{y) 

Gi (pi(y)) - Pg (pi(y)) tr (V5(y) - id). 



(60) 



Combining - (I^U)) . we conclude that 

G{poiy))dy- [ G{pi{y))dy > 



Pcipiiy)) tr (yS{y)- id) dy 
Pg {pi{y)) div {S{y) - y) dy. 



Next, assume that Pg{pi) G W^'°°{^) and pi > a.e. Since Pg > 0, we can approximate 
Pg{pi) by nonnegative functions in C^{M'^). We use Proposition 14.11 - (iv), to obtain 
that 



PG{pi{y))div{S{y)-y)dy > I {V[PG{pi{y))], S{y) - y) dy (61) 

(V[G'(pi(y))],5(y)-y)pi(y)dy. 
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We combine (|54p and (|61l) to conclude (|5()|) □ 

The next theorem extends the energy inequality (|56jl to general cost functions c. 
Theorem 2.8 (Energy inequality for general cost functions). 

Let po, Pi £ T^ai^) be such that pi > a.e., and c : R'^ — > [0, oo) be strictly convex, 
of class and satisfy c(0) = and (H(J^. Let G : [0, oo) R be differentiable on 
{0,oo), such that G(0) = 0, (0,cxd) 3 x \—>- x'^ G{x~'^) be convex and nonincreasing, 
V{G'{pi)) £ L°°(17), and Pcipi) G W^'°°{n). Denote by S, the c-optimal map that 
pushes pi forward to pQ. Then, 

Gipo{y))dy- [ G{pM)dy> [ {V[G' {pM)l S{y) - y) pi{y) dy. (62) 
Jn Jn 

Proof : Let {ck)k be a sequence of regular cost functions satisfying (|36|) . By Proposition 
12.71 we have that 

G{po{y))dy- [ G{pi{y))dy> [ {V {G' {p,iy))) , Sk{y) - y)piiy) dy, (63) 
Jn Jn 

for all k £ iV, where denotes the Cfc-optimal map that pushes pi forward to po. We 
let A: go to cxD in (|63j) . and we use that V{G'{pi)) £ L°°(ri), and {Sk)k converges to 
S in Ll^{n,R'^) ([2], Lemma 2.2.2), to conclude (jHH); here L'j,^{n,R'^) denotes the set 
of measurable functions : 17 — > R'^ whose square are summable with respect to the 
measure pi := pi{y) dy, i.e. | ip{y) P/Oi(y) dy < oo. □ 

2.4 Approximate solutions to the paraboHc equation 

Throughout this section, we assume that po + £ L°°(r2). For fixed h > and i £ iV, 
we denote by the minimizer of 

(P,^) : inf {hW^{pl„p)+E,ip) : p £ P„(0)}, (64) 

where p^ := po (see Proposition 12. 3|) . We define the approximate solution p'* to ©, as 

r po{x) if t = 
p''it,x):=l (65) 
[ pUx) if t£ [ti^uti], 

where ti = ih, i £ iV. The next proposition shows that 



dp^ 



div{/ Vc* [V (f'(p^))] } +A(/i), 



dt 

in a weak sense. We show in the next section, that 

||A(/l)||(^2,oo(n)). =0(/l<'')) 

where e{q) := min(l, q — 1). 



21 



Proposition 2.9 Assume that F : [0, oo) R is strictly convex and satisfies F G 
C^((0,oo)), and c : R'^ [0, oo) is strictly convex, of class C^, and satisfies (Hd^ - 
(H(i^. Then 



T 



Jn 



{PQ - p^) d'li dx dt + / I {p'^Vc* V(F'(/)) ,VC)dxdt (66) 



Jo. 



1 

^ [o,T]xn 



T/h 

D'i{t,x)\ / 



where, ^ : R x n ^ R is such that ^{t, .) G C^{n) for t G R, and spt(,{.,x) C [-T,T] 
for X £ Q, and for some T > 0. Here, 



^{t + h,x)-^{t,x) 
h 



and 7j is the Ch-optimal measure in T{pf_^, p^) . 

Proof : Without loss of generality, we assume that f G iV. Because of (gSI), we have 
that 



Ai{t,x)dx 



< B' 



for t G (0, T), where 



Af{t,x) :-- 



h 



at,x) + lp^{x)Vc^ viF'ip^ix)]] ,Vat,x)), 



and 



i.f:=- sup 



D''^it,x) I \x-y\^d-iKx,y). 
Jnxfi 



[o,T]xn 

We integrate the above inequality over t G (0, T), to obtain that 



dt A'^{t,x)dx < h^Bt 

i=l Jii-i j=l 



(67) 



The right hand side of 1)6 7p gives that 

D^at,x)\Y[ \x-y\^dj^{x,y), 

i=l •'f^xf^ 



T/h 



sup 



T/h 



(68) 



i=l - [0,T]xQ 

while, on the left hand side, we have that 



1=1 



T/h 



Y / / A^{t,x)dxdt = Y 



1=1 



+ 



pKx) - P'l-iix) 



ti—l ''f^ 

T 



Cit,x)dxdt (69) 



Jn 



p'^Vc* 



V(F'(/) 



) dxdt. 
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By a direct computation, the first term on the right hand side of ()69|) gives that 



/^^^^^^^^e(^,x)dxd^ = l£ ly{t,x)at,x)dxdt 

1 ^ r 

7-E / / p\r-h,x)ar,x)dxdT 



_2 li-i 

1 



, , / Poix)£,it,x)dxdt. 
We use the substitution t = t + h in the above expression, to obtain that 



Iff If 

- / p^it,x)^{t,x)dxdt - - p^{t,x)^{t + h,x)dxdt 

Jo Jn " Jo 



/ Po(2;) x) dx dt 


T f T rT 

h 



[ [ p''{t,x)d'^^{t,x)dxdt + \ f 
Jo Jn " Jt 



p^{t,x)C{t + h,x) 

T-h 



h r- 

I Po{x) ^{t,x) dtdx. 
Jn 



Noting that 



T f f P0ix)Cit,x)dtdx= [ [ poix)d^C{t,x)dxdt, 

Jo Jn Jo Jn 

and ^{t + h) = ioi t £ (T — h, T), we deduce that 



Jn 



E/ / /^^ ^ ;'-^^-^^ e(x,t)dxdt (70) 



T 

Pq{x) - p^{t, x) ] d^^it, x) dx dt. 



10 Jn 

We combine dHZI) - dZOl) to conclude □ 

3 Existence and uniqueness of solutions 

Below, we study the limit of H66|) . as h goes to 0. The first three sections deal with 
the limits of the three terms of inequality (|66|) . and the last section proves the existence 
theorem to ©, when pQ is bounded below and above. Here and after, p^ is defined as 
in (jHSI). 
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3.1 Second moments of the optimal measures 

We show that 

T/h 

|x-ypd7,"(x,y) = 0(/i<«)), (71) 

where e{q) := mm(l,g — 1), denotes the c/i-optimal measure m r{p^_^, p'l), and p-* 
is the unique minimizer of (|64j) . The first step toward proving (|71|) is the next lemma, 
which states that Yli^i ^ciPi-i^ Pi) is bounded, uniformly in h. 

Lemma 3.1 Assume that F : [0, cxd) M and c : R'^ [0, 00) are strictly convex, and 
c satisfies (Hd^. Then 

p^hW^ipl^p'l) < E^{po) -\n\F . (72) 

Proof: Let T > 0, be such that ^ G N. Since c(0) = 0, Proposition 12.31 gives that 

hW^{pl,,p1)<Ei{pl,)-E,{p^), 
for i G IV. We sum both sides of the above inequality over i, to obtain that 

T/h 

Y,hW^{pU,p^)<EM- / F(p^/,(x))dx. 
1=1 •'^ 

We apply Jensen's inequality to the integral term above, and we let T go to 00, to 
conclude (O □ 

Proposition 3.2 Assume that F ■.[Q,oo) ^ R and c : R'^ ^ [0, 00) are strictly convex, 
and c satisfies c(0) = and c{z) > (3\ z\^ , for some (3 > and q > 1. Then, for T > 
and h G (0, 1) such that ^ G W, 

T/h 

Y \x-y\^d^^{x,y)<Min,T,F,po,q,f5)h'^'i\ (73) 

where e{q) := min(l,g — 1). 

Proof . Since c{z) > P\z\'^, we have that 

hi 

\x-y\idjH^,y) < -tW^{pI,,p1), (74) 

for i G iV. We distinguish two cases, based on the values of q. 
Case 1: I < q < 2. 

Because of ()74() . we have, for i G IV, that 

/ |x-ypd7,^(x,y) < sup jx-yl^^--?) / \x - y\^ dj^{x,y) 
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where, diamO denotes the diameter of fi. We sum both sides of the above inequality 
over i, and we use (f7^. to conclude that 



T/h 

E 

i=l 



x-y\^ dj^x, y) < Min, F, po,q, f3) h^'K 



Case 2: q > 2. 

Because of Jensen's inequality and (fTl]) . we have that 



nxn 



X - yl"^ d'yi{x,y) < 



\x-y |M7f(x,y) 



2/q 



< 



2/q 



We sum both sides of the above inequality over i, and we use Holder's inequality on the 
right hand side term, to obtain that 

2/q 

U2 /rr\ i-^^ 

X -y\^d-/i{x,y) < 



T/h 

E 



nxn 



01 



T/h 
i=l 



T/h 
i=l 



2/q 



(75) 



□ 



We combine (|7^ and (f7^. to conclude that 

T/h 

\x-y? d-i^{x, y) < M{n, T, F, po, g, (5) h 

3.2 Strong convergence of the approximate solutions 

We prove that {p^)h is compact in L^{^}t), for < T < oo. The main ingredient in 
the proof is the energy inequality (|62() . It allows us to obtain a uniform bound in h, of 
the L'^ - norm of V , which leads to the compactness of {p^)h in L^{^t)- We 

first show that {p^)h converges weakly in L^{^t) for a subsequence. We introduce the 
following constant needed in the next lemma: 



M{n, T, F, PO, q, a) := M(a, q) Ei{po) - \n\F 



+ aT\^ \ II Po ll^o 



where M{a, q) is a constant which depends only on a and q. 

Lemma 3.3 Assume that c : [0, oo) is strictly convex, of class and satisfies 

(H(^), and F : [0,oo) ^ R is strictly convex, of class ((0, oo)). // po G Va{^) n 
L°°{n), then, 

II llL°°((0,oo);L°°(f7)) < II Po \\L°°{n)- (76) 

Therefore, there exists p : [0, oo) ^ M and a subsequence of (p^O/iiO; which converges 
to p, weakly in V'{U,t), for < T < oo. 

In addition, if ^ ^ L°°{Q), c satisfies (H(^, and F satisfies F{0) = and (HI^, then 



V(F'{p^' 



< M{n,T,F,po,q,a). 



(77) 
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Proof : Because of the second bound in (|26|) - a consequence of the maximum prin- 
ciple of Proposition 12.21 -. we have that < \\ pQ \\L°°(^n) for i G N, which reads as 
II p^{t) \\L°°{n) < II Po \\L°°{n) for t G [0, oo). We take the supremum of the previous in- 
equahty over t £ (0, oo), to deduce ((7^ . 

Due to (f76|) . we have that {p^)h is precompact in L^{^It), for < T < cxd. We use the 
standard diagonal argument, to conclude that {p^)hio converges weakly to some function 
p : [0, oo) xVL ^ M\n L^{ilT)-, for a subsequence. 

Because of Proposition EHl (EEl), and the fact that V {P{p^)) = p^V {F'{p^)), we have 
that P{p^) G VFi'°°(rj) and V (F'(pf)) G L°°(fi), for i G iV. Then, we choose G := F 
in the energy inequality (jMJi and we use (jl^ . to obtain that 



V(nP.')),Vc^ V(F'(pf)) )pf</F(pf_i) 



np^)- 



We sum both sides of the subsequent inequality over i, and we use Jensen's inequality, 
to deduce that 



V (F'(p'^)) ,Vc* V 



/< / F(/,o)-|f^|F(^ 



1 



(78) 



Because of H139|) of Proposition 14.21 and the fact that c[z) < a (| z 1"^ + 1), we have that 

(z,Vc*(z)) > c*{z) > M{a,q)\zf - a, 
and then, ((75)) implies that 

1 





V (f'(/)) 








Jn 



F{po)-\n\F( — ]+a 



P 



(79) 



We combine (|75)) and (f7^ . to obtain that 





V (f'ip")) 











^F(po)-|J^|f(^) +aT|l]|||/5o 



L°°(n)- 



We divide both sides of the above inequality by M{a, q) to conclude ((77|) 



□ 



Lemma 3.4 (Space-compactness) 

Assume that c : R'^ [0, oo) is strictly convex, of class and satisfies c(0) = 
and (H(][^, and F : [0, oo) M is strictly convex, of class C^((0,oo)), and satisfies 
F{0) = and (H^. If po G Va{^) is such that Po + j^^ L'^i^), then, for r] and 
< T < oo, 



/ p''{t,x + rje)-p\t,x) <Min,T,F,po,a,q)\vl 



where e is a unit vector of K^^ Q.^"^^ := {a; G : dist{x,dO,) > | ?7 1 }, and 0,> 
(0,r) X 



(80) 



M ._ 

T ■- 
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Proof: Since /Oq + ^ G L°°(J7), (PH]) implies that {p^)h is bounded below and above. 



Then, we use that F ^ C ((0, cxd)), to obtain that 



L9*(Ct) 



1 



VlF'ip- 



< M{n,po,F) 



ViF'ip^) 



(81) 



We combine (|77)) and (|5T|) . to conclude that (y P^)j^ is bounded in {^t)- As a con- 
sequence, we have that {p^)h is bounded in T^^''^ {Q.t)- Approximating p^ by C'^{Q.t)- 
functions, and using the mean-value theorem, and the fact that {'Vp^)h is bounded in 
L"^ (f^r), we have that 

[ \p\t,x + r]e)-p''{t,x)f <Min,T,F,po,a,q)\rif. (82) 
We combine ()82() and Holder's inequality, to conclude that 



('?) 



p^{t,x + 'ne) - p''{t,x) 



(')) 



p^{t,x + r]e) -p^(t,x) 



< M(l^,T,F,po,a,g)|r/ 



□ 



Now, we focus on the time-compactness of {p^)h on VLt-, < T < oo. The following 
constant will be needed in the next lemma: 

^(17,r,F,po,9,«,/3) 

II n 11^/9* 

II P'Ollioo(Q) 



L°°(n) 



M(g, a, /3) ( EM - \n\F [ — ] + aT\n\\\ po ||loo(o) 



where M(g, a, /3) is a constant which only depends on g, a and /3. 

Lemma 3.5 Assume that c : R'^ [0, cxd) is strictly convex, of class and satisfies 
c(0) = and (H(^, and F : [0,oo) M is strictly convex, of class ((0,cxd)), and 
satisfies F(0) = and (H^. If po £ Va{^) is such that po + £ L'^^i^), then, for 
r > and < T < cxd, 

" F' (p\t + r, x)) - F' (p\t, x)) 1 [ p\t + T,x)- p\t, x) 



<Min,T,F,po,q,a,l3)T. 



(83) 



Proof : Without loss of generality, we assume that ^ £ N and T = Nh, for some 
N £ W. For simplicity, we set 



Lih,T) :-- 



F' [p\t + T, x)) - F' [p\t, X)) ] [p\t + T,x)- p\t, x) 
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and 

It is straightforward to check that 

T/h 

L{h,T) = ^hJ{i,h,N). 
1=1 



(84) 



Since (Wj})^/'^ does not satisfy the triangle inequahty, we introduce the g-Wasserstein 
metric := (W^'')^/'', defined by 



y-sHy) 



(85) 



where denotes the | ^ | ''-optimal map that pushes p^_^j^ forward to p^. Then, setting 
^In P' [p^+n) - [p'i)^ we obtain that 



J{i,h,N) 



<7v(y)-<7V pf+Niy)dy 



Since Po + ^ G L^^m, F e ((0, cx))), and p^V {F'{p1)) = V {P{p^)) G L-(0) 
(see Proposition EH - (i)), ^ gives that <^^^ G W^'°° (Q) . So, approximating v^^^ 
by C°°(r2)-functions, and using that (5^)^/9^^^ = p^ and the mean- value theorem, we 
rewrite J{i, h, N) as follows: 

J{i, h,N)= V^l^ ((1 - t)y + tS.^y)) , y - S'^{y) ) p^M dtdy. 

We combine Holder's inequality and (|85|) . to deduce that 

j(i,/i,iv)</i4(p^ptiv) / y'|v(/.tiv((i-i)y + <(y)) |%tiv(y)dtdy 

^(86) 

But, observe that p^, p^j^^ < || po llL°°(r2) because of Proposition 12.21 and | V(/?^^ \'^* G 
L°°(r2). So, we approximate | V(/9^^ |''* by nonnegative functions in C^{R'^), and we 
use in Proposition 12.41 to deduce that 

V^^.':^ ((1 - t)y + tp^(y)) P ptiv(y) dy < || Po ||L-(n) | V(^,';^(y) [ dy. (87) 

We combine (jHHl) and (|HZ1), to have that 



L{h, r) < II PO Ili(!;^) ^ 4 (p?, Ptiv) II V^f,^ IL, 



T//i 



'(C)- 



i=l 



And since is a metric, the triangle inequality gives that 

j_ N T/h 

L{h,T) < llpollll(f^) h^^Y,\\^'Pi,N\\L'>*{n)dg [p'l+k-i^Pi+k 

k=l i=l 
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Then, we apply Holder's inequality to the interior sum, to deduce that 



L{h,T) < II po 



/T/h 



Because of © and (|77|), ( h^/'^* V {F'{p^)) 



Li*{n) 



N 

E 

k=l 



i=l. 



T/h 
1=1 



and 



1/9 



h 



belong to lq*{i}). Then, we combine Holder's 



inequality, Minkowski's inequality, ()26() and (|77() . to have that 



T/h 

E" 




^( 




^1 


EH 

Vi=i 



g 

LI* (Q) 



LI* (n) 



+ 



V(F'(pf 



L9* (Q) 



rN, 



1 

Li*{Q) 



1/q* 



1 

Li* {Q) 



1 



^ [M{n,T,F,po,q,a) 



1 1/9* 



1/9 



On the other hand, since c{z) > P \ z\'^, we have that (dg)'^ < ^W^, and then. 



N 

E 

k=l 



T/h 



i=l 



1/9 



N 



1 



/3' 
T/h 

i=l 



1/9 



We use (|72|) and the above inequality, to deduce that 



N 

E 

k=l 



T/h 



1/9 



1 

< - 



Ei(/)o)-|l^|F 



17 



1/9 



Nh-^''i. (90) 



We combine (|55)) - (jHOJ, and we use that r = A^/i, to conclude that 

L(/i,t) <M(17,r,F,po,g,a,/3)r. □ 
Lemma 3.6 (Time- compactness) 

Assume that the assumptions of Lemma \3.4\ hold. If po G Va{^) is such that po + ^ £ 
L°°{0,), then, for < T < oo, and small r > 0, 

/ p\t + T,x)- p^{t,x) <M{R,n,T,F,po,a,q,(3)V^ + TAiV^), 
where A is such that limT-j^g ^(\/^) = 0- 
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Proof : Let R > 0, and for fixed h, T and r, define 



F' {p\t)) , + 



+ 



F' (/(t + r)) 



11(0) 



Because of dTT)), and the fact that F £ C^{{0,oo)), we have that (0,r) 3 

t ^ Ah,T{t) belongs to ((0,r)). Hence 



Er \ < 



M{n,T,F,po,q,a,f3) 



R 



We combine (|7^ and (|MT) . to have that 



/ / p^(t + r,x)-/(t,x) <2||/)o||L-(n)|^^||£^i?,| < 



M(0,T,F,po,g,a,/3) 



(91) 



(92) 



On the other hand, if t E := (0, T) \ ^^j, setting /(t) := ui and p''(t + r) : = 
M2, we clearly have that || lij ||L<j(n) < II -^'(""j) llvFi''?*(f2) < for i = 1,2, and 
[F'{u2) — F'{ui)] [u2 — ui] < Rt. Then, Proposition 12.51 gives that 



[ f \p^{t + T,x)- p^{t,x)\< f A{Rt)<TA{Rt) 



(93) 



where A{Rt) := Ar{Rt) is defined as in ProDosition l2.51 We combine (|^^ - and 
we choose R = to conclude the proof □ 

Having proved the space-compactness and time-compactness of {p^)h7 we are now 
ready to show that {p^)h converges strongly to p in L^{^t), < T < oo, for a subse- 
quence, where p is defined as in Lemma 13.31 

Proposition 3.7 Assume that c : R'^ [0,oo) is strictly convex, of class and 
satisfies c(0) = and (H(l[^, and F : [0, oo) R is strictly convex, of class ((0, oo)), 
and satisfies F(0) = and (H^. If pQ G Va{^) is such that /)o + ^ G i°°(^^), then, for 
< T < oo, there is a subsequence of {p^)hio which converges strongly to p in U{VIt), 
for 1 < r < oo, where p is defined as in Lemma Vj/A 

Proof: Fix 5 > O, and define fi^'* as in Lemma 13.41 Because of H76() . we have that 
{p^)h is bounded in (^^^^^ ■ Furthermore, for e > 0, and small r > and 77 G (0, 5), 
we have that fi^^ C Jl^^ C and then. Lemma 13.41 and Lemma 13.61 give that 

/ I /(t, x + rie)- p^{t, x)\<e, and / | /(t + r, x) - /(t, x) | < e. 
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uniformly in h. We deduce that, (yo'*)h is precompact in ^^t^^ (^^^ EP' Theorem 

2.21). We observe that hm^^o I ^ \ ^^^^ I = 0) and then, we use the diagonal argument, 
to obtain that {p^)h converges strongly to p in L^{^t), for a subsequence. And since 
{p'')h is bounded in L°°(Qt) (see we conclude that it converges to p in 

for 1 < r < oo (up to a subsequence) □ 

3.3 Weak convergence of the nonlinear terms 

We use the energy inequality ^ to show that (div{p^Vc* [V {F' (p^))]}) ,^ converges 
weakly to divjpVc* [V (F'(^))]} in ily, for a subsequence. Throughout this section, 
{p'^)h denotes the (non-relabeled) subsequence of {p'^)h which converges to p in L'^{^t)^ 
for 1 < r < oo, as in Proposition 13.71 and 



h ._ 



V(F'{p^) 



The next lemma shows that {cr^)h is bounded in L'^{Q,oo), and (V (F'(p^)))^ converges 
weakly to V {F'[p)) in {^t) for a subsequence.. 

Lemma 3.8 Assume that c : R'^ — > [0, oo) is strictly convex, of class and satisfies 
c(0) = and c{z) > (3 \ z\'' for some (3 > and q > 1, and F : [0, oo) R is strictly 
convex, of class C"^ ((0, oo)). If Pq € Va{^) is such that po + ^ S L°°(r2), then 



\rr^\\'^ < 



/3 



po 



Eiipo 



1 



(94) 



(i). Therefore, there is a subsequence of ('^'')^jo' '^^'^^^ converges weakly to a function 
a in L'^i^r), for <T < oo. 

In addition, if c satisfies (II(][^, and F satisfies -F(O) = and (HI^, then 



(a), there is a subsequence of |v (yF'{p^)) | , which converges weakly to V {F'{p)) 
in L"* ((0, T) X n), forO<T <oo. 

Proof : By we have that 

sHy)-y 



V(F'(pf(y)) 



(95) 



for i € iV, where denotes the c/i-optimal map that pushes forward to Pi_i. We 
use (|26p and (|95|) . to deduce that 



h \\q 



i=l 



< 



Y^h lJvc*[v[F'ip1iy)))] 'dy = Y,h 

sHy)-y " 



i=l 



sHy)-y 



dy 



, . i=i ■'^ 



J_ 




po 


L°°(Q) 



Pi{y)^y- 
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Since c{z) > P \ z\'^, we obtain that 



IIl9(Qoo) - 



j_ 

PC) 



(96) 



i=l 



We combine (|7^ and to conclude (i) is a direct consequence of (j^lj. 

Now, fix < T < oo. By Proposition 13.71 {p^)h converges strongly to p, in L^{Qt), 
and by ((7^ and the fact that F' is continuous on (0, oo), (^F'{p^)^^ is bounded in 
L°°(Ot)- We deduce that [F'{p^))^ converges weakly to F'{p) in L'>* (Qt)- And, since 
I V {F'{p^)) 1^ is bounded in L'^* (Qt) (because of JSI and ((13)), we conclude (ii) □ 

The next lemma extends the energy inequality 1)62(1 to the time-space domain (0, oo) x $7. 
Lemma 3.9 (Energy inequality in time-space) 

Assume that c : [0,oo) is strictly convex, of class C^, and satisfies c(0) = 

and (H(^, and F : [0,oo) R is strictly convex, of class C^((0, oo)), and satisfies 
F(0) = 0, F G C2 ((0,oo)) and (H^. If po G Va{^) is such that Po + & L'^i^), 
and t u{t) is a nonnegative function in C'^{M), then 



n 
1 

< - 
- h 



o^V {f'{p^) 



,Vc' 



V [F'ip'^ 



F(po(x))n(t) + 



u{t) 
F{p'^)d^u{t), 



where 



d^uit) :-- 



u{t + h) - u{t) 
h ■ 



Proof : Let T be such that ^ G IV, and assume that spt u C [— T, T] . We choose G := F 
in the energy inequality and we use to obtain that 

F {pHv)) - F {ptAy)) 



< 



h 
V 



dy 



F' [p1{y)) ] , Vc* [V [F'{p^{y))) ] ) pf(y)dy. 



for all i G iV. Since u > 0, we deduce that 



T/h 

E 

i=l 



/■ F(pf(y))-F(pti(y)) 



u{t) 



< 



/(V(F'(/)),Vc^ ViF'ip^)) )uit). 
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By direct computations, the left hand side of the above inequality gives that 

[ F {p'iiy)) - F {pl,{y)) 



T/h 

E 

i=l 



ti-l •'^2 



■u{t) 



(97) 
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F(/(t,y))^(t) - F{po{y))u{t) 
F[p\t-h))u{t). 



IT 

f-T 



h Jn 



We use the substitution t = t — h in the last integral, and the fact that u{t + h) = for 
t € {T - h,T), to have that 

fZMM^pftM,^,^ (98) 

F{p\t,y)) dMt)-l [ F{po{y))u{t). 

We combine (^7)) and (P5|) . and we let T go to oo, to complete the proof □ 

Theorem 3.10 Assume that c : M'^ — > [0, oo) is strictly convex, of class C^, and sat- 
isfies c(0) = and (H(J^, and F : [0, oo) R is strictly convex, of class ((0, oo)), 
and satisfies F(0) = 0, and (HI{^ - (H^^. If po G P„(S7) is such that Po + j^^ L'^i^), 
and t H-> u{t) is a nonnegative function in C'^{]R), then 

hm / ( A", V (f'(/)) ) u{t) = j {pa,V{F'ip)) > u{t), (99) 

where p and a are defined in Lemma and Lemma 

Therefore, (^div{p^a^)) ^ converges weakly to div{pa) for a subsequence in [C^{M x R'^)'j , 
and 

div{pa) = div{pVc* [V {F'{p))]) . (100) 

Proof : Let T > be such that sptu C [—T,T], and assume that p(t) = po, for t < 0. 
Denote by {p^)h the subsequence of {p^)h, such that 

(i) . {p^)hio converges to p a.e., 

(ii) . { V (F'(/)) }hio converges weakly to V{F'{p)) in L'^^Qt), and 

(iii) . {cr^ = Vc* [V (F'(p^))] }hio converges weakly to a, in L1{VLt), 
as in Proposition Em and Lemma 13.81 We first observe that 

lim/ {a\p^V{F'{p)))u{t)= [ {a, pV {F' (p))) uit), (101) 

and 

hm / (p^Vc'^ [V (F'(p)) ] , V (f'(p'^)) - V {F'{p)))u{t) = 0. (102) 

Indeed, since (/);, is bounded in L°°(Qr) (see and V (F'(p)) G L^^VIt), (i) 

and the dominated convergence theorem imply that {p^V {F' {p))}hiQ converges to 
pS/ {F'{p)) in L"? {^t)- Then, we use (iii) and the fact that u G C^{R), to conclude 

((TnT|i . 
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Because of Proposition 14.21 the convexity of c, and the fact that c{z) > P \ z\'^, we have 
that 

|Vc*(^)|«<£lZ^ = ^{{z,Vc*{z))-c%z)) 

< ^{z,Vc^z)) <M{P,q)\zf. 



We deduce that 



Vc*[V(F'(p))] <M{(3,q) V(F'(p)) 



(103) 



which shows that Vc* [V {F'{p))] G L^(r2j'). Then, we use (i) and the dominated con- 
vergence theorem, to have that { p^Vc* [V {F'{p)) ] }hio converges to pVc* [V {F'{p)) ] 
in L'?(Qt)- We conclude p()2|) . because of (ii). 
The proof of H99|l fohows directly from the following three claims: 

Claim 1. 

/ {pa,V {F'ip)) > uit) < hm inf / ( V ) u{t). 

Proof : Because c* is convex, and u and p^ are nonnegative, we have that 

\ p^ (Vc^ [V (f'(p^))] - Vc^ [V [F\p))\ , V {f\p^)\ - V [F\p)) ) > 0, 
and then, 

hminf/ (a^p'^V(F'(p)))n(^) 

<liminf (A',v(F'(p'^)))^x(t) 

+ limsup / (p'^Vc* [V ] , V - V {F\p'^)])u{t) . (104) 

We combine (|101|) - ()104() . to conclude Claim 1. 
Claim 2. 

limsup / //cj\V ('r(p''))\ u(t) 



Proof : First, we observe that 



< / [poF'{po)-F^F'{po))] u{0) 
Jn 

+ [ [pit,x)F'{p{t,x))-F^{F'{p{t,x)))]u'{t). 



hm / F{p')d^u{t)= [ F{p)u'it) 



(105) 
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Indeed, it is clear that 

F{p^)d^u{t)-F{p)u'{t) 

[ \Fip^)-F{p)\\u'{t)\+ [ \Fip^)\\d^u{t)-u'it)\. (106) 



< 



Because of and the continuity of F, we have that [F{p^)^^^ is bounded in L°°{^t)- 
We let go to in ()1U6() . and we use (i), the fact that u G Cl{R), and the Lebesgue 
dominated convergence theorem, to conclude H1U5|) . 
Lemma EiH gives that 

limsup / {p^a^\V (f\p^))) u{t) 

<liminfi/ [ F{po)u{t)+ limsup [ F{p^)d^u{t), 
MO h Jo Jq Jq^ 

and by ()1U5() and the continuity of u, we deduce that 

limsup / (p^a\V (F'ip''))) u{t)< [ F{po)uiO)+ [ F {p{t, x)) u' (t). (107) 

Since F £ ((0, oo)) is strictly convex and satisfies F{0) = and (Hl[T]), we have that 

F* {F'{a)) =aF'{a)-F{a), Va > 0. (108) 
We substitute H108|) into 1)107(1 for a = p(t,x) and a = po{x), to conclude Claim 2. 



Claim 3. 



[poF'{po)-F* {F'ipo))] n(0) 

+ / [pit,x)F'{p{t,x))-F^{F'{p{t,x)))] u'{t) 
It 

{pa,V{F'ip)))u{t). 



< 



Proof: Set ^(t,x) := F' {p{t,x)) u{t) for {t,x) e R x n. Because of (i), (ii), (HU, and 
the fact that F e ((0,oo)), we have that F'{p) G L°°(Jlr) and V {F'{p)) G LI^Qt)- 
We approximate F'{p) by C°°(r2T)-functions in W^''^ {^t), and we use with the 
backward derivative d^^ S,{t,x) := , and Proposition 13.21 to obtain that 

/ {p,-p>-)d;'^i+ I {a\p'^y{F'{p)))u{t)=m^^\ 

where t{q) = min(l,g — 1). We let /i go to in the subsequent equality, and we use 
(|101|1 . to conclude that 

hm/ {p,-ph)dr^C+ [ {a,pV{F'{p)))u{t)=0. (109) 
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Since sptu C [— T, T], we have that 



poix)S,{t,x) 



and then, 



lim / podr''^ 



We combine (fTUUI) . ((TTU)) and (i), to have that 



PoF\po)uiO). 



(110) 



(a,pV(F'(p)))n(t)=lim / x)5,A(t, + / poF' {po)u{0) . (Ill) 

By direct computations, we obtain that 

p{t,x)d^''at,x) = p{t,x)F' {p{t,x))d^^u{t) 



h 



p{t, x) u{t -h)[F' {p{t, x)) - F' {p{t - h, x)) ] . 



Since F ^ ((0, oo)) is strictly convex, and satisfies F[0) = and (Hl[T]), we have that 

{F'{b) - F'{a)) b > F* {F'{h)) - F* (F'(a)) , V a, 6 > 0, 
and then, we deduce that 



pit,x)d;-''at,x) > p{t,x)F'{p{t,x))dtK{t) 

+ 1 u(t -h)[F^ {F' {p{t, x))) - F* [F' {p{t -Kx)))]. 

We integrate both sides of the subsequent inequahty over $1^, and we use that n = on 
(T — h,T) for h smaU enough, and p{t,x) = po{x) for t £ (— /i,0), to obtain that 



p{t,x)dt''i{t,x) > / [p{t,x)F'{p{t,x))-F^{F'{p{t,x)))]dr''u{t) 



I [\{t-h) I F* (F'(po)). 



We let /i go to in the above inequality, to deduce that 



lim 



p{t,x)d^''i{t,x) > / [p{t,x)F'{p{t,x))-F*{F'{p{t,x)))]u'{t) 



F^ (F'(po))n(0). 



(112) 



We combine ()111() and (|112() . to conclude Claim 3. 

In the end, we show that a = Vc* [V {F'{p)) ], which combined with Lemma 13.81 com- 
pletes the proof of Theorem 13.101 Indeed, let e > 0, ip £ C°^{Q), and set LOe{t,x) := 
F' {p{t, x)) — €ijj{x). It is clear that Vu^ £ L"? (Qt), and 



Vc*(Vw,) < M{P,q)\VuJ,\^ , 
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as in the proof of (|1U2|1 . We deduce that Vc* (Vcj^) G L'^{Qt)- We use that c* is convex, 
and and u are nonnegative, to have that 



(Vc* [ V (f'{p^)) 1 - Vc* (Vw,) , V - Vw,) u{t) > 



We let /i go to in the above inequahty, to obtain that 

hmsup / {p^ a^, V (F'{p^)))u{t) -limM [ {a^ , p^Vw,) u{t) 

-hminf y (/Vc*(Vu;,),V (^F'(/)) - Vw,) u(t) > 0. (113) 

As in the proof of H1U1|) and 1)102^ . we have that 



hminf / {a^ , p^Vwe) u{t) = {a, pV-w^) u{t), (114) 



and 



■ Vtfe) n(t) 



(115) 



limM [ (p^Vc*(Vw,),V (f'(p^ 
hlo JnT ^ \ yr 

= [ {pVc*iVw,),V{F'{p)) -Vw,)uit). 
We combine ^ and (|TT^ - (fTT^ . to have that 

/ {pa - pVc*(Vwe),V (F'(p)) - Vwe) u{t) > 0. 
We divide the subsequent inequahty by e, and we let e go to 0, to obtain that 

{pa-pVc* [V (F' (/>))] ,VV(x)n(t)) >0. 
Choosing —ip in place of ip, we get that 

/ {pa - pVc* [V (F'ip)) ] , V^{x) u{t)) = 0. 

JQ.T 

And since ^ and n > are arbitrary test functions, we deduce H1UU|1 . This completes 
the proof of the theorem □ 

3.4 Existence and uniqueness of solutions 

Here, we state and prove the theorem of existence and uniqueness for 0. 

Theorem 3.11 (Case V = Q). 

Assume that c : — > [0, cxd) is strictly convex, of class C^, and satisfies c(0) = 
and (Hdl^, and F : [0,cxd) ^ R is strictly convex, of class C^((0, oo)), and satisfies 
F{0) = and (H^ - (H^. If po G Va{^) is such that Po + ^ G L'^i^), and V = 0, 
then, 0) has a unique weak solution p : [0, oo) x $7 ^ [0,oo), in the sense that 
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(i). P+-p€ L°°((0,oo);L°°(17)), V G L''*(17t) forO<T <oo, and 

(ii). forCeC^{Rx R'^), 

L {-P^ + ip"^^* ['^{F'ip))]^^0} = I^Poi^)^(^^^)d^- (116) 

Proof : Proposition 13 . 71 gives that {p^)h converges to p a.e. for a subsequence, and since 
p'^ > for all h, we deduce that p > 0. We combine H^IEh and Proposition 13.71 to have 
that /0 + i G L~ ((0,oo);L°°(l^)). We use that V(F'(p)) G L'^* (Qt) (LemmaEHI to 
conclude (i). 

Recall that Xui^ gives that Vc* [V (F'(/j))] G L'?(Ot) for < T < oo, and ^ and 
Proposition O imply that p G L°°(17t). We deduce that pVc* [V (-F'(/o))] G L5(0t)- 



Now, fix < T < oo, and let ^ G 



be such that spt^(.,2;) C [-T,T] for 



x G ri. Because of Proposition 12.91 and Proposition 13.21 we have that 



lim 



(Po 



a,'^e+(p"Vc'^[v(F'(p'^))],VO}=0. 



(117) 



Lemma [3.31 gives that {p^)h converges weakly to p in L^(Ot) for a subsequence, and 
then, we have that 



lim 



{po 



(118) 



From Theorem Eini we have (div{/Vc* (V {F'{p^)))})^ converges weakly to 
div{pVc* (V {F'{p)))} in [C^{R x J?"')] , for a subsequence, then we deduce that 

hm/ {p'^Vc^[V{F'{p))],VO= [ (pVc^[V(F'(p))],VO. 
HO Jq Jq^ 



(119) 



We combine (|117() - (|119() . and we use the fact that spt^(.,x) C [— T, T], to conclude 

jnni). 

Here, we prove uniqueness of solutions to © when ^ G ((0, T) x fi), for < T < cxd. 

Using the arguments in jEI, it is easy to extend the proof to the general case. In fact, 

assumption ^ imposed in |3] would not be required here. The convexity of c*, that is, 

{Vc*{zi) — Vc*{z2), z\ — Z2) > 0, for zi, 2:2 G 1?*^, suffices to extend the proof. 

Let T > 0, and assume that pi and p2 are weak solutions of © with the same initial data, 

such that N < pj < M a.e., and ^ G L^^^t), j = 1,2. Since V {F'{pj)) G Li*{nT), 

and 

Vc* [V (F'ip,)) ] P < M{P,q) I V (F'(p,)) 
we have that V c* [V {F' (pj)) ] G L9(17t). For 5 > 0, we define 

f^T B {t, x) ^ ^5 {t, x) := ifs [F' {pi{t, x)) - F' {p2{t, x))) , 

where 



if T < 
5 if 0<T < 6 



1 if T > (5. 
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Using a smooth approximation of S^g as a test function in the differential equations 
satisfied by pi and p2, and passing to the hmit, we have that 

/ i5dt{pi-p2) = - f (piVc*[V(F'(pi))] -/>2Vc^[V(F'(p2))],VC5), 
which reads as 

/ £.5dt[pi-p2) 

-] - '^2)(Vc* [V {F'{p2)) ] ,V{F'{p^) - F'{p2))), 

where n>fy := n[0 < F'{pi) - F'{p2) < S]. Because c* is convex, the first term on 
the right hand side of the above equality is nonpositive. And since F S ((0,oo)) is 
strictly convex and satisfies (HF^, and N < pi, p2 < M a.e., we have, a.e. on ^2^^ , 
that 



\Pi- P2\ = 

We deduce that 

isdtipi - P2 



[{F*y o F'] (pi) - [{F^y o F'{p2)] < 6 sup (Frir). 

re[F'(Af),F'(M)] 



< sup (Frir) 

T(^[F'{N),F'{M)] 



(1.2) 



{Vc*[V{F'{p2))],V{F'ip,)-F'{p2))) 



(1 2) I 

We let 5 go to in the subsequent inequality, and we use that tps — > -^[o,oo)) I ^t's I ~^ 
and [F'(pi) - F'(/)2) > 0] = - p2 > 0], to have that 



dt[{pi-p2)^] <0, 



Pi(r)-p2(T)]+< / [pi(o)-p2(o)]+ = o, 



which reads as 



for < T < oo. Interchanging pi and p2 in the above argument, we conclude that 
Pi= P2 □ 



Theorem 3.12 (General case). 

Assume that V : Q ^ [0, oo) is convex, of class C^, and c : — > [0, cxd) is strictly 
convex, of class , and satisfies c(0) = and (H(J^. Assume that F : [0, oo) — > R 
is strictly convex, of class C^((0,oo)), and satisfies F{0) = and (Hl^ - (Hl^. 
If po e Va{^) is such that po + ^ £ F°°{^), then, ^ has a unique weak solution 
p : [0, oo) X ^ [0, oo), in the sense that 
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(i). P+-p€ L°°((0,oo);L°°(17)), V(F'(p)) G L''*(17t) forO<T <oo, and 
(ii). for e G C2(J? X 1?'='), 

/ {"'"M^^'"^"' [V(F'(p) + y)],VC>}= / po(x)^(0,x)rfx. (120) 

Proof : The proof of the uniqueness of solutions is similar to that of Theorem 13.111 
Here, we only prove existence of solutions to (jH)). Let ^ € C^{R x R'^) be such that 
spt^(-,x) C [— T, r] for X G and for some < T < oo. Following the arguments in 
the previous sections, where the minimization problem (P) (|16|) is replaced by 



(P^) : ini{hW!,^{po^p)+E{p) : peVain)}, 



and 



E{p) ■.= E,{p)+ [ pVdx, 
Jn 

we have, as in Proposition 12.91 that 



{po-p'')d^Cdxdt+ / (/Vc* V(F'{p^) + V) ,VOdxdt 

T/h 



1 

< sup 

^ [0,T]xQ 



D^i{t,x)\Y,l \x-y\''d^Hx.y). (121) 

4 = 1 •'f^Xf^ 



and, as in Proposition 13.21 that 



T/h 

E 

i=l 



y |2 d7f (x, y) < M(17, T, F, po, 9, /3) /i^^^^ • 



We let /i go to in ()121() . and we use (|122() . to deduce that 

hm / {(po - p") d^i + { p'^Vc* [ V (f'(p^) + f) j , VO) = 0. 
The following claim suffices to conclude Theorem 13.121 



Claim. For < T < oo, the estimates 



\P llL°°({0,oo);L°°(f7)) < II PO ||L°°(n)) 



(122) 



(123) 



(124) 



P 



<Min,T,F,po,V,q,a), 



V [F'ip^) 

and the energy inequality in time-space 

(p^V (f'{p^) + v),Vc*\v (f'{p^) + V) ) u{t) 



< 



[F{po)+poV]uit)+ Fip^)+p'^V d^u{t) 



(125) 



(126) 
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hold, for nonnegative functions u in C^{R). 

Indeed, because of (|124|1 . there exists p : [0, oo) x Q, ^ [0, oo), such that 
(iii). {p^)h converges to p, weakly in L}{VIt), for a subsequence. 
As a consequence, 

hm / (po-/)a,H= / (po-p)§. (127) 

Using H124|) and (|125() . we deduce the space-compactness and the time-compactness of 

{p'^)h in L^{^t), as in the case where V = 0. Hence, 

{iv). {p^)h converges strongly to p in L}{Q.t), for a subsequence. 

Then, we use (if), (|126|) . and we follow the lines of the proof of Theorem I3.in( where 
we use F'{p^) + y in place of F'{p^), and F{p^) + p^V in place of F{p^), to conclude 
that 

{v). (div{p^Vc* [V (F'(/) + V)]})^ converges weakly to div{pVc* [V {F'{p) + V)]}, 

in [Cl{R X -R'')]', for a subsequence. 
Hence, 

hm/ (p'^Vc^[v(F'(/) + y)l,VO= / (pVc^[V(F'(p) + y)],VO. (128) 
We combine (fT^ . (ITTTl) . and dUHl), to conclude (fnUj) . 

As in Theorem 13. Ill (i) follows directly from (jl24|) . (|125|) . and the maximum/minimum 
principle of Proposition 12.21 for \/V ^ 0. 

Proof of the Claim : ()124() is a direct consequence of the maximum principle of Propo- 
sition O for VV / 0. 

As in the case F = 0, we have, because of Proposition 12.61 and the maximum/minimum 
principle of Proposition O that P{p^) E W^^°°{Vl), and V (F'(pf )) G L°°(S1). Then, 
choosing G := -F in Theorem 12. 8| the (internal) energy inequality (|(i2j) read as 

/ F{pl,)- [ F{(/:)> [ {v(F'{p'^)),sHy)-y)pHy)dy, 

Ja Jn Jn ^ 

where is the c/^-optimal map that pushes p^ forward to Pi_i- We use that {Sj^)^p^ = 
Pi_i, and V G C^{^) is convex, to deduce the potential energy inequality 



pIiV- j/iy^ I {yv,sHy)-y)pHy)<dy- 



We add both of the subsequent inequalities, and we use the Euler-Lagrange equation of 
(P^), that is. 



h 



V{F'{p\{y)) + V{y)) , for a.e. y G J^, (129) 



(where is the c/i-optimal map that pushes p\ forward to po)) to deduce the free energy 
inequality 

E{pl,) - E{p^) >hJ^{V [f'{p^) + y) , Vc^ [ V {f'{p^) + v)]) pt (130) 
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for i G N. We sum (|13Uj) over i, and we use that V and p^^^ are nonnegative, and 
Jensen's inequality, to have that 



(v(F'(/) + y),Vc^ v(F'(/) + y) )p'^<Eipo)-\n\F( — ]. 



We conclude, as in the proof of ((77|) . that 



V(F'(p^) + F) ' <M(17,r,F,po,g,a)- 



On the other hand, because of H124|) and the fact that V £ C^{Q), we have that 



Li*{nT) 



< II Po 



|l/<3* 



\L°°{n)- 



(131) 



(132) 



We combine (fnT|) and (HSU, to conclude (fT^ . 

The proof of (|126jl follows the lines of the proof of Lemma 13.91 where we use the free 
energy inequality H13U|) in place of the internal energy inequality (|62() □ 

Remark 3.13 (Existence of solutions to ^ for a wider class of pq). 
Here, we show how to extend our existence theorem to a wider class of initial 

probability densities po : po £ L°°(r2) and ^ L°°(r2), or po G LP{Q), p > q. For 
simplicity, we assume that V = 0. 

Case 1: po G L°°(rj), and ^ L°°(f7). 

Let {po^s)5 be a sequence in Va{^), such that 

m < Po,s < II Po ||L°°(n) a.e., where < ■i]s < 5 
Ei{po,5) < Eiipo) 

K Po,5 Po in L^{^), as 5 i 
(see [2], Proposition 1.4.2), and define the approximate solution p^ to by 

/9o,5 if t = 
pI, if tG ((i-l)/i,i/i], 
where p^^ is the unique minimizer of 

(P,,5): inf{Wi^(pti^„p)+i?,(p): p G Pa(0)}. 
Since /9o,5 + G ^°°(f^)) we have as before, that 



(133) 



P\--- 



II P5 llL°°((0,oo);L°°(n)) < II ^0,5 ||L°°(n)) 

. Kp\ ^{E'{p\)) ' <M{n,T,F,po,s,q,a), 



(134) 
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and 



(135) 



where e{q) := min(l, q-1), C^{B x W^), and 
M{n, T, F, po,5, q, a) := M{a, q) ( Ei{po,s) -\n\F 



+ aT\n \ II po,5 ||l° 



We introduce a convex function H : [0, cxd) R, such that, 

(HHl) : (7i([0,oo))nC2((0,oo)) and = xi/'?*_p"(2;), Vx > 0. 

Combining ([HH]) . (fTT?H) . and (HHU), we have that 

I ^5 llL°°{{0,oo);L°°(r2)) < II PO IIl°°(Q)) 



Li*(nT) 



< M{n,T,F,po,q,a). 



(136) 



We deduce that there exists p : [0, oo) x $7 ^ [0, cxd), such that {p^)h,5 converges to p in 
Li(Qr) for a subsequence, and {a^ := (p^)i/'?Vc* [ V (F'(p^)) ] j?,,^ converges weakly to 
pi/^Vc* [V(F'(/3))] in L'^iyiT) for a subsequence, as (/i, <5) goes to (0,0). Then, we let 
{h,5) go to (0,0) in ()135() . to conclude that p G -L°° ((0, oo); L°°(r2)) is a weak solution 
of ©, as in Theorem 13.111 but here, we do not require that ^ G ((0,cx));L°°(Q)) □ 

Case 2: po £ ^'^{Vl), p> q, and Ei{pQ) < co. 

Using Corollary 1.4.3 |2|, we approximate po by a sequence {po,s)s in 'Pa(^), such that 
% < /Oo,5 < a.e., where < % < (5 and £5 > | 

^i(po,5) < ^i(po) 

. Po,s Po in LP (17) as (5 i 0, 

and we define /O^ as in case 1. Since po L°°(i7), we cannot obtain (|136j) from (|134|) . 
as in case 1. Here, we take advantage of the fact that po £ LP{Q) and Ei{pQ) < 00, as 
follows: 

(i). We choose G{x) := x > 0, in the (internal) energy inequality and we 



observe that (V//'-^, Vc* [V {F'{p))]) > for p G Vai^), to have that 



(137) 



II Ps l|L°°({0,oo);LP(n)) < II Po \\LP{n)- 

As a consequence, there exists a function p : [0, 00) x ^ [0, 00) such that {p^)h,5 
converges weakly to p in L'p{Q.t), < T < 00. 

(ii). Next, we choose G := F in to control the spatial derivatives of p^, as 
P^\^{F'{pt 

<M{a,q) (E,{po)-\n\F 



1 



Po \\Lp{n) 



(138) 
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We combine (|137|) . l|138|) . and we use (HI^, to deduce, as in the previous sections, that 
{Ps)h,5 converges strongly to p in L^{Qt), for a subsequence. We conclude, as in case 1, 
that /) is a weak solution of © in the sense that, p £ L°° ((0, oo); LP(r2)), 
pVc* [V {F'{p))] G L^i^T) for < T < oo, and ((TTHl) holds □ 

Few examples of energy density functions satisfying (H_b[T]) - (HFl^I) and (HI^ are 

= E"=i AiFi{x), where Fi{x) £ {xln(a;), with m > 1 or max 1 - < 

m < 1, and Ai > 0. For examples, for the fast diffusion equation ^ = Ap"*, this corre- 
sponds to the range 1 — ^<m<lifd>2, and ^ < m < 1 if d < 2; for ^ = App"^, we 
have n > '^^(^Si^ d > p, and n > p(^p__i-j d < p, and in particular, for the p-Laplacian 
^ = Ap/9, we require < p < d or p > max (^d, ^-'^^^ • 

4 Appendix 

In ProDosition l4.11 we collect results of previous authors used in this work, and in Propo- 
sition we establish intermediate results needed in the previous sections. Proposition 
14.11 is due to Cordero [S] and Otto ^Sj . For its proof, we refer to these references. A 
sketch of proof of this proposition can also be found in [2j, sections 5.1 and 5.2. 

Proposition 4.1 Let po,pi £ Va{^), and assume that c : R'^ — > [0, oo) is strictly 
convex, and satisfies c,c* £ C'^{R'^). Denote by S, the c-optimal map that pushes pi 
forward to po, and define the interpolant map St, and the interpolant measure pi-t, by 

St := {I - t)id + tS and pi-t '■= {St)#pi, 

for t £ [0, 1] . Then, 

(i) . St is injective, and pi-t is absolutely continuous with respect to Lebesgue. 

Moreover, there exists a subset K of Q, of full measure for pi := pi{y)dy, such 
that, for y £ K and t £ [0, 1], 

(ii) . VS{y) is diagonalizable with positive eigenvalues. 
(Hi). The pointwise Jacobian det{SJS) satisfies 

/ Pi {y) = Pi-t {St{y)) det [{l-t)id+ tVS{y) ] , 
where pi-t is the density function of pi^f 

In addition, if pi > a.e., then 
(iv). the pointwise divergence div{S) is integrable on Vl, and 

div {S{y) - y) i{y) dy < - [ { - y, VO dy, 
n Jn 

forC>0 in C^{R'^). 
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The following estimates will be needed in the previous sections. 

Proposition 4.2 Assume that c : R'^ — > [0, oo) is strictly convex, of class , and 
satisfies c(0) = and (H(^. Then 

{z,Vc''{z)) >c*{z)>0, yz e R'^. (139) 

In addition, if c{z) > (3\z\'^ for some (3 > and q > i, then 

{z,Vc*{z)) < M{f3,q)\zf , (140) 

where M{p,q) is a constant which only depends on [i and q. 

Proof : Since c is strictly convex, differentiable and satisfies (HCEJ, we have that c* G 
C^{R'^) is convex. Then, 

(z, Vc* (z)) = c*(z) + c (Vc*(z)) > c*(z). (141) 

Because c(0) = and minimizes c, we have that c*(0) = and minimizes c*. We 
conclude that c*{z) > 0, which proves (|139)) . 

Now, assume that c{z) > j3\z\'^. Since c* G C^{R'^) is convex and nonnegative, we have 
that 

( z, Vc*{z) ) < c*{2z) - c*{z) < c*(2z). (142) 
Moreover, because c{z) > (3\z\'^, we have that 

^{2z) <M{l3,q)\zf . (143) 

We combine (fTl^ and (fTlH|) . to conclude (|Tin|l □ 
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